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We solve the Dirac equation approximately for the attractive scalar S{r) and repulsive vec- 
' tor V{r) Hulthen potentials including a Coulomb-like tensor potential with arbitrary spin-orbit 

coupling quantum number k. In the framework of the spin and pseudospin symmetric concept, 
Q^, we obtain the analytic energy spectrum and the corresponding two-component upper- and lower- 

' spinors of the two Dirac particles by means of the Nikiforov-Uvarov method in closed form. The 

limit of zero tensor coupling and the non-relativistic solution are obtained. The energy spectrum 
for various levels is presented for several k values under the condition of exact spin symmetry in 
I the presence or absence of tensor coupling. 
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I. INTRODUCTION 



Within the framework of the Dirac equation the spin symmetry arises if the magnitude 
of the attractive Lorentz scalar potential S{r) and the time-component repulsive vector 
potential are nearly equal, S{r) ~ V{r) in nuclei {i.e., when the difference potential A(r) = 
V{r) — S{r) = Cs = constant). However, the pseudospin symmetry occurs when S{r) ~ 
—V{r) arc nearly equal (i.e., when the sum potential S(r) = V{r) + S{r) = Cps = constant) 
[1-3]. The bound states of nucleons seem to be sensitive to some mixtures of these potentials. 
The cases A(r) = and E(r) = correspond to SU{2) symmetries of the Dirac Hamiltonian 
[3]. The spin symmetry is relevant for mesons [4]. The pseudospin symmetry concept has 
been applied to many systems in nuclear physics and related areas [2-7]. Further, it is used 
to explain features of deformed nuclei [8], the super-deformation [9] and to establish an 
effective nuclear shell-model scheme [5,6,10]. The pseudospin symmetry appeared in nuclear 
physics refers to a quasi-degeneracy of the single-nucleon doublets and can be characterized 
with the non-relativistic quantum numbers {n,l, j — I + 1/2) and {n — 1,1 + 2, j — 1 + 3/2), 
where n, I and j are the single-nucleon radial, orbital and total angular momentum quantum 
numbers for a single particle, respectively [5,6]. The total angular momentum is given as 
j = I + I, where I = I + 1 is a. pseudo-angular momentum and s = 1/2 is a pseudospin 
angular momentum. In real nuclei, the pseudospin symmetry is only an approximation and 
the quality of approximation depends on the pseudo-centrifugal potential and pseudospin 
orbital potential [11]. The Dirac Hamiltonian with vector and scalar potentials quadratic in 
space coordinates has been studied [12]. It was shown that the the Dirac equation can be 
solved exactly for the cases A(r) = and E(r) = 0. In addition, the linear tensor potential 
and mixture of quadratic scalar and vector potentials are studied [13]. It is shown that 
a linear tensor potential with quadratic A(r) or S(r) generates a harmonic-oscillator-like 
second order differential equation which can be solved analytically. Recently, Akcay [14,15] 
has shown that the Dirac equation for scalar and vector quadratic potentials including 
the Coulomb-like tensor potential with the spin and pseudospin symmetries can be solved 
exactly. These results in Dirac equation with quadratic potential plus a centrifugal-like 
potential can also be solved analytically. Tensor coupling potentials are added as spin-orbit 
coupling terms to the Dirac Hamiltonian by making the substitution — imulS.rU (r) 

[16,17]. Tensor couplings and exactly solvable tensor potential have been used to investigate 
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nuclear properties [18-20] and have also some physical applications [21,22]. 

In the past years, there has been much interest in the solution of the relativistic Dirac 
and Klein-Gordon equations [1,12,23-26]. For instance, some authors have solved these 
equations for several physical potentials, such as the Woods-Saxon potential [23-26], the 
Morse potential [27], the Hulthen potential [28], the Eckart potential [29-31], the Poschl- 
Teller potential [32,33] and the Scarf- type potential [34], etc. 

Recently, many works have been done to solve the Dirac equation so to obtain the energy 
equation and the two-component spinor wave functions. Jia et al [35] employed an improved 
approximation scheme to deal with the new centrifugal spin-orbit term k{k-\-1) in the 
second order differential equation that results from the Dirac equation and to solve it for 
the generalized Poschl- Teller potential for arbitrary spin-orbit quantum number k. Zhang et 
al [36] solved the Dirac equation with equal Scarf-type scalar and vector potentials by the 
method of the supersymmetric (SUSY) quantum mechanics, shape invariance approach and 
by alternative methods. Zou et al [37] solved the Dirac equation with equal Eckart scalar and 
vector potentials in terms of SUSY quantum mechanical method, shape invariance approach 
and function analysis method. Wei and Dong [38] obtained approximately the analytical 
bound state solutions of the Dirac equation with the Manning- Rosen for arbitrary k. Thylwe 
[39] presented the approach inspired by amplitude-phase method in analyzing the radial 
Dirac equation to calculate phase shifts by including the spin- and pseudo-spin symmetries 
of relativistic spectra. Alhaidari [40] solved Dirac equation by separation of variables in 
spherical coordinates for a large class of non-central electromagnetic potentials. Berkdemir 
and Sever [41] investigated systematically the pseudospin symmetric solution of the Dirac 
equation for spin 1/2 particles moving within the Kratzer potential connected with an angle- 
dependent potential. Recently, we have also solved the spin and pseudospin symmetric Dirac 
equation with arbitrary spin-orbit centrifugal term for generalized Woods-Saxon potential 
[42] and Rosen-Morse potential [43] by means of the Nikiforov-Uvarov (NU) method. 

In this paper, it is worth to investigate the solution of the Dirac equation for scalar and 
vector Hulthen potential for A(r) or S(r) together with Coulomb shape tensor coupling 
potential which can be solved analytically by using an improved approximation scheme 
introduced in Refs. [44,45] to deal with the resulting centrifugal and pseudo-centrifugal 
terms k{k±1) r~^. The Coulomb-like tensor potential preserves the form of the Hulthen 
potential but generates a new spin-orbit centrifugal terms A (A ± 1) r~^, where A is a new 
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spin-orbit quantum number. This provides a possibility for generating a different form of 
spin-orbit coupling which might have some physical applications. 

The Hulthen potential, widely used for the description of the nucleon-heavy nucleus 
interactions, takes the following form (see [44] and the references therein): 

Vnir) = -^^r^, = r\ V, = Ze\ (1) 

where Vq is the potential depth, 5 is the screening range parameter and Tq represents the 
spatial range. If the potential is used for atoms, then Vq = Z5 (in the relativistic units 
h = c = e = 1), where Z is identified as the atomic number. The Hulthen potential behaves 
like the Coulomb potential near the origin (i.e., r — )> or r <^ ro) Vc{r) = —Ze^/r , but 
decreases exponentially in the asymptotic region when r ^ 0, so its capacity for bound states 
is smaller than the Coulomb potential. This potential has been applied to a number of areas 
such as nuclear and particle physics [46-48], atomic physics [49,50], molecular physics [51,52] 
and chemical physics [53], etc. 

In the presence of the spin and pseudospin symmetry, we investigate the bound state 
energy eigenvalues and corresponding upper and lower spinor wave functions for arbitrary 
spin-orbit k quantum number in the framework of the NU method [54-56]. We also show 
that the spin and pseudospin symmetric Dirac solutions can be reduced to the S{r) — 
V{r) and S{r) — —V{r) in the cases of exact spin symmetry hmitation A(r) = and 
pseudospin symmetry limitation E(r) = 0, respectively. Furthermore, the solutions obtained 
for the Dirac equation can be easily reduced to the Schrodinger solutions when a parametric 
transformation is applied. 

In what follows, we first review the NU method and present a parametric generalization 
in Sect. 2. In the presence of spin and pseudo-spin symmetries, we obtain the bound state 
solutions of the Dirac equation with scalar and vector Hulthen potentials including the 
Coulomb-like tensor interaction, the limit of zero tensor coupling and the non-relativistic 
limits by applying a suitable transformation in Sect. 3. The relevant concluding remarks 
are given in Sect. 4. 
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II. THE NIKIFOROV-UVAROV METHOD 



The NU method [54] is briefly outhned here. It is based on solving the second-order 
differential equation of hypergeometric-type by means of special orthogonal functions: 

r{r) + ^V^'W + ^Mr) = 0, (2) 

where a{r) and a(r) are polynomials at most of second- degree, 'r(r) is a first-degree polyno- 
mial and %l){r) is function of the hypergeometric type. In order to find a particular solution 
for Eq. (2), we choose 'il){r) as follows: 

^(r) = (j){r)y{r), (3) 

which leads to a new second-order hypergeometric-type equation of the form 

y"{r) + A{r)y'{r) + B{r)y{r) = 0, (4) 

with A{r) and B{r) are taken to be T{r)/a{r) and a'{r)/a'^{r), respectively, where a(r) = 
A(T(r), A is a constant and r(r) is a polynomial of degree at most one. If we impose that 
4>{r) satisfies the following logarithmic equation 

0(r) (7(r)' ^ ^ 

then we obtain 

r(r) = T(r) + 2n{r), T'{r) < 0, (6) 

where 7r(r) is a polynomial of order at most one. Making use of (p" {r) / (pir) = (0'(r)/0(r))'-|- 
{(j)' {r) / (j){r)Y , we can reduce Eq. (4) into another hypergeometric-type: 

a{r)y"{r) + T{r)y'{r) + Xy{r) = 0, (7) 

and the quadratic equation for the polynomial 7r(r), 

7r^(r) + [f (r) - a'{r)] n{r) + a{r) - ka{r) = 0, (8) 

where 

A; = A-7r'(r). (9) 
The solution of the above quadratic equation for 7r(r) yields 

7r(r) = ^ [a'{r) - f (r)] ± ^1 [^'(r) -T{r)]' -a{r) + ka{r), (10) 



and the weight function can be obtained via 

[(7(r)p(r)]' = r{r)p{r), 



(11) 



where the prime denotes the differentiation with respect to r. The expression under the 
square root sign in Eq. (10) can be arranged as a polynomial of second order where its 
discriminant is zero. Hence, an equation for k is being obtained. After solving such an 
equation, the k values are determined through the NU method. One, however, is looking 
for a family of solutions corresponding to 

X^Xn^-nT'{r)-^n{n-l)a"{r), n = 0,l,2,--.. (12) 

The y{r) — yn{r) which is a polynomial of degree n can be expressed in terms of the 
Rodrigues relation: 

yr^^'^ = ];(7)d^n[-''ir)p(r)h (13) 

where Bn is the normalization constant and the weight function p(r) in (13) is the solution 
of the differential equation (11). 

In addition, an eigenvalue solution through the NU method can be set up from the 
relationship between A and A„ by means of Eqs. (9) and (12). 

For a more simple application of the method, we develop a parametric generalization of 
the NU method valid for any central and non-central exponential-type potential by making 
change of the independent variables z — z{r). Thus, we obtain another generalized hyper- 
geometric equation 

[z (1 - C3z)f w"{z) + [z{l- C3z) (ci - C2z)] w' {z) + {-q2Z^ + qiz - qo) w{z) = 0, (14) 

when compared with Eq. (1) yields 

t{z) = ci- C2Z, a{z) = z{l- c^z) , a{z) = -q2Z^ + qiz - qo, (15) 

where the parameters Cj and qj {j — 0, 1, 2) are to be determined during the solution 
procedure. Thus, by following the method, we can also obtain all the analytic polynomials 
and their relevant constants necessary for the solution of a radial wave equation. Hence, 
these analytical expressions are displayed in Appendix A. 
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III. DIRAC EQUATION WITH A TENSOR COUPLING 



In spherical coordinates, the Dirac equation for fermionic massive spin-1/2 particles mov- 
ing in attractive scalar S{r), repulsive vector V{r) and a tensor U{r) potentials reads as (in 
relativistic units h — c — 1) [13,57]: 

[a • p + ^ (M + S{r)) + V{r) - ipa-rU{r) - E] V'n«(r) = 0, Vn«(r) = V'(r, 9, 0), (16) 

where E is the relativistic energy of the system, M is the mass of the fermionic particle, 
p = — i V is the three-dimensional momentum operator, and a and (3 are 4x4 Dirac matrices, 
which have the following forms, respectively. 









(l \ 


a — 








[a 1 




[O-I 



(17) 

where I denotes the 2x2 identity matrix and a are three- vector Pauli spin matrices 

<7i= I ° M , a2= I ° ~M , a3= I ^ ° I . (18) 





For a particle in a spherical (central) field, the total angular momentum operator J and the 
spin-orbit matrix operator K = — /3 (o" ■ L + I) commute with the Dirac Hamiltonian, where 
L is the orbital angular momentum operator. For a given total angular momentum j, the 
eigenvalues of K are k = —{j + 1/2) for aligned spin (si/2, P3/2, ^tc.) and k = j + 1/2 for 
unaligned spin (pi/2, c?3/2, etc.). The spinor wavefunctions can be classified according to the 
radial quantum number n and the spin-orbit quantum number k and can be written using 
the Pauli-Dirac representation: 

ipnK[r,0,(l)) = - \ , (19) 

where Fn^ir) and GnK{r) are the radial wave functions of the upper- and lower-spinor com- 
ponents, respectively, Yjj^{9,(f)) and Y-^^{9^(j)) arc the spherical harmonic functions coupled 
to the total angular momentum j and it's projection m on the z axis. The orbital and 
pseudo-orbital angular momentum quantum numbers for spin symmetry I and pseudospin 
symmetry I refer to the upper- and lower-components, respectively. For a given spin-orbit 
quantum number k = ±1, ±2, • • • , the orbital angular momentum and pseudo-orbital an- 
gular momentum are given hy I — \n-\- 1/2| — 1/2 and I — \k — 1/2| — 1/2, respectively. 



The quasi-degenerate doublet structure can be expressed in terms of a pseudo-spin angular 
momentum 's — 1/2 and pseudo-orbital angular momentum I which is defined as I — I -\- 1 
for aligned spin j = 1 — 1/2 and I — I — 1 ior unaligned spin j = 1 + 1/2. For example, 

(3si/2, 2^3/2) and (3^1/2,2^3/2) can be considered as pseudospin doublets. 
Substituting Eq. (19) into Eq. (16) and using the following relations [57] 

(cr- A) (cr-B) = A -B + ifj- (A X B), (20a) 

((7 • P) = (7 • r ( r • P + ) , (20b) 



and properties 

{a-h)Y]^{e,<l>) = {K-l)Yf^{e,<l>), (21a) 
{a ■ L) YlJe, 0) = - (« + 1) YIJ0, 0), (21b) 
(a-f)l^L(^,0) = -l^U^,0), (21c) 
(a-r)Y^e,<l>)^-Yj^(e,<l>), (21d) 
we obtain the following two radial coupled Dirac equations for the spinor components: 

{/l + ^- U{r)^ Fn.{r) = (M^ + K« - A(r)) Gn.{r), (22a) 

{jr~1^ t^(^)) GUr) = {M' - E^^ + S(r)) F,,(r), (22b) 
where A(r) = V{r) — S{r) and S(r) = V{r) + S{r) are the difference and sum potentials, 
respectively. By ehminating GnK{f) in Eq. (22a) and -F„«;(r) in Eq. (22b), we get two second- 
order non-linear differential equations for the upper and lower radial spinor components, 
respectively 

^ + 1) ^ _ _ ^ . ^,^^1 ^^^^^^ 



dr"^ r dr M -|- — A(r) \dr r 

= (M + B„, - A(r)) (M - B„, + E(r)) F„.(r), (23) 



(ir^ r M — £'„k^ -|- E(r) \dr r 

= (M + E,, - A(r)) (M - E,, + S(r)) ^^.(r), (24) 

where k (k -|- 1) = l{l + 1) and k,{k, — 1) — l{l + 1). The radial wave functions are required 
to satisfy the necessary boundary conditions, that is, F„k(0) = Gn/t(0) = and ^^^(r) = 
GnK{f) — >■ at infinity. 
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At this stage, we take the E(r) or A(r) the form of Hulthen potential (1) and the tensor 
potential in the form of the Coulomb-like interaction. Equations (23) and (24) can be solved 
exactly for k — 0,-1 and k — 0,1, respectively, because of the spin-orbit centrifugal term. 



A. Spin symmetric bound state solution 

In this part we are taking the E(r) as the Hulthen potential and A(r) = Cs — constant 

E(r) = V„(r) = -^^7^, S„ = V„, (26) 

In fact, when the limit of Tq becomes infinity, then lim T,{r) — oo, it tells us that the Dirac 
particle could not be trapped by the Hulthen potential, which does not have any bound 
state under the condition of spin symmetry. The Coulomb-like potential [58] for the tensor 
due to a charge Z^e interacting with a charge Zab, distributed uniformly over a sphere of 
radius Rc, is added, 

Ucouiir) = --, r > R,, (26) 

where Rc — 7.78 fm is the Coulomb radius, and Za and Za denote the charges of the 
projectile a and the target nuclei A, respectively. Under this symmetry, substituting Eqs. 

(25) and (26) into Eq. (23), the equation obtained for the upper radial spinor Fn^{r) becomes 



' ' ■F„,(r)=0, (27) 



(f 1 , (3 £nk 



(^j,2 J.2 (e^Ao _ 1) j-^ j-^ 
where 

l^VkiVk-i), Vk^ K- + H + 1, (28a) 

13 = 4 {En. + M-Cs)Vo> 0, (28b) 

^nk — '"o (K« - M) {Enn + M-Cs)< jdo, M > En., (28c) 

with K — I and ac = — (/ -|- 1) for k > and ac < 0, respectively. Also, the quantum condition 
in (27) is obtained from the finiteness of the solution at the origin point 

FnM = 0, (29) 

and at infinity 

limF„,(r) = 0, (30) 
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The spin symmetric energy eigenvalues depend on n and k, i.e., Enn — E{n, k {k + 1)). 
In Eq. (25), the choice of E(r) = 2V{r) — >■ V{r) as mentioned in Ref. [12] allows us to 
reduce the resulting relativistic solutions into their non-relativistic limits under appropriate 
transformations. 

Equation (27) can not be solved exactly only for the case of k = 1 due to the spin-orbit 
centrifugal term 7r~^ which is expanded in terms of singular functions of e~'^^^° compatible 
with the solvability of the problem for r <^ tq. Now, since the orbital term is too singular, 
the validity of such approximation is limited only to very few of the lowest energy states. 
Therefore, to go to higher energy states one may attempt to solve the relativistic version 
of the problem in [44] like the Dirac equation (27). Therefore, wc apply the approximation 
scheme derived in [44] for the centrifugal term which is valid for small screening parameter 
S values (i.e., for large ro values or r <^ ro). It can be casted in the form [44]: 



+ 



er/ro _ 1 ^gr/r-o _ 1)^ 



(31) 



where the dimensionless constant do — 1/12 is exact as reported in many recent works (cf. 
e.g., [35,59-62]). Obviously, the above approximation to the centrifugal (pseudo-centrifugal) 
term turns to when the parameter ro goes to infinity (small screening parameter 5) as 



lim 

ro— >-oo 




sr/ro 



which shows that the usual approximation is the limit of our approximation (cf. e.g., [44] 

and the references therein). In terms of the new variable z(r) = e"^''''", which maps the 
interval r e (0, oo) into z G (0, 1), and using the approximation in Eq. (31), then Eq. (27) 
transforms into rational functions as 



{1- z) d 
z{l — z) dz z^ 



where 



qo = 'ydo - £nk = A^^, 



(32) 

(33a) 
(33b) 
(33c) 
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and XnK must be a positive real parameter for the presence of bound states. We are looking 
for solutions in the form 

= z^{l - zYfiz), e > 0, ^ > 1, (34) 

where ^ and <; are real positive parameters and the boundary conditions in Eqs. (29) and 
(30) are also satisfied. Then, by substituting (34) into (32) the following hypergeometric 
differential equation for f{z) is obtained 

z{i - ^)^/(^) + [2e + 1 - (2e +2^+1) z] j-j{z) - [(e + 0' - Q2] m = o, (35) 

Thus, the wave functions satisfying Eqs. (27), (30) and (34) are given by 

FnAz) - M^«(l - z^'^z^"'^ {-71, H + 2A„« + 2rjk + 1, 2A„« +l;z), (36) 

where 



1/1 , \ 1 + H +1, K> 0, 

= o + \/ 7 + ?o + ?2 - gi = « + + 1 ^ 77^ = < (37) 



-l + H, K<0. 



2 V 4 

For the bound-state problem, the solutions (36) fulfill the boundary condition (29) when 

AnK + 77fe- -ri, n = 0,1,2,- •• , (38) 

where n is the number of the nodes of the radial wave functions. This also results in the 
following energy spectrum formula 

P-jn + t^ + H + lf 
^"'^ - 2{n + K + H + l) ■ ^^^^ 

The last formula gives an equation for the energy levels. When we insert A^^ and (3 into the 
above equation, we find the following transcendental energy spectrum formula 

f/o (/V + H) (/V + H + 1) 



{En. - M) {Enn + M-Cs) 





2 



rp {En. + M-Cs)Vo {n + K + H + 1) 
{n + K + H + 1) To 



n = 0, 1, 2, • • • ,n„,ax = r-oVVo {Enn + M - Cs) - k - H - 1, (40) 

where rimax being the largest integer which is less than Tq-^/Vq {En. + M — Cg) — k — H — 1, 
Enn > Cg — M. In the hmit of zero tensor couplings {H — 0), the spin aligned states {k < 0) 



11 



demand that n + 1 ^ —k. However, such a condition is no longer exist for the spin unahgned 
states > 0). 

On the other hand, we may also apply the NU method, this can be done when Eq. (32) 
is compared with Eq. (14), the following polynomials are found 

t{z) = 1-2;, cr{z) = z{l - z), a{z) = q2z'^ - q^z + go- (41) 

Furthermore, we can follow the parametric generalization model of the NU method given in 
Appendix A to obtain the specific values for the constants q {i — 1, 2, • • • ,13), the results 
are listed in Table 1 for the potential model under consideration. When the relations (Al- 

A4) of Appendix A together with the values of constants given in Table 1 are applied, the 
key polynomials take the following particular forms [63]: 



(42) 



and 



A; = /3 - 7 - 2A„, [Vk-\]: (43) 



for discrete bound state solutions. According to the NU method, we further obtain 

Tiz) = 1-2 



r'(z) = -2(^% + A„, + <0, (44) 

with prime denotes the derivative with respect to z. Referring to Table 1 and applying the 
relation A5 of Appendix A, we obtain energy equation formula which is identical to Eq. 
(40). 

In the limit of zero tensor coupling and under the usual approximation, H — do — Cg — 
K — I and Vo — Ze^5, Eq. (40) becomes 

which is completely identical to Eq. (68) of Ref. [64] for the solution of the Klein-Gordon 
equation with equally mixed scalar and vector Hulthen potentials, i.e., Vq — Sq (Eq = 0). 
Furthermore, the last equation is also identical to Eq. (25) upon inserting q — 1^ Vq — So 
and S — S± — 1 + 1,-1 in Eqs. (8) and (9) of Ref. [65]. Therefore, real solutions are possible 
only for \Eni\ < M (i.e., bound states). 
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We now look at some special cases and relationships between our results and some other 
existing results in literature. For exact spin symmetry, Vq = -S'o case or Cs — and applying 
the following appropriate transformations Eni + M ~ 2/x, E^i — M ~ E^i and k — I, we 
obtain the following non-relativistic energy spectrum for the Hulthen potential including the 
Coulomb-like interaction, 



Eni = ^i{H + l) {H + l + l)do 



fiZc 



_S{H + n + l + l) 2 
leading to the following energy spectrum formula in the limit of zero tensor coupling [44] 



(45) 



E.^ 



6' 



nl 



1(1 + 1) do 



fiZ( 



{n + I + 1) 



2M \ ' ' ^ S{n + l + l) 2 
n = 0, 1, 2, • • • , rimax = ro^/2jJ,Vo - I - I, 



(46) 



which is identical to Eq. (34) in Ref. [44] for the solution of the Schrodinger equation with 
the Hulthen potential. The largest integer nmax < '"ov^2/IVo — / — 1 which is identical to Eqs. 
(12) and (13) given in Ref. [66] after choosing A = 2firlVo, a = 1 and a = for the Hulthen 
potential case, respectively. In fact, when the limit of S becomes zero, then 

lim Enl = — — TT — — ^ — oi 

r-o^oo 2M{n + l + lf 

which is a spectrum resulting from a Coulombic field. 

We apply now the relations (A6-A10) of Appendix A to calculate the corresponding wave 
functions as [63] 

p(^) =^2A„«(^_^)2^«^ (47) 

0(^) = ^^'-(1 -;2)''«, (48) 

where 77^ > 1 and A„„ must be a positive real parameter for the presence of bound states. 
Hence, we find 

Vniz) = p(^^^A)^i -2z), ze [0, 1], (49) 

with Pn''^\x) is the Jacobi polynomial, where z/>— l,/x>— 1 and defined for the argument 
X e [—1,1]. Using the relation FnK{z) — (f){z)yn{z), we get the radial upper-spinor wave 
functions as 

n!r(2A„, + 1) 



T(n + 2A„« + 1) 
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(50) 



Furthermore, we give the relation hnking the hypergeometric function and the Jacobi poly- 
nomials (see formula 8.962.1) in [67] 

1-x 



Fi -n, n + u + iJ, + 1,1' + 1; 



to rewrite the radial wave functions as 



r(n + i/ + l) " ^ ^' 



(51) 



Fn.{r) = Nn.e-^^-''''^'- (1 - e-^/'")"^ ^F^ (-n, n + 2A„, + 2r;, + 1, 2A,, + 1; e"''/'") . 

(52) 

where the normalization constant is calculated in Appendix B in closed form. 

Before presenting the corresponding lower-component Gni^{r), let us recall a recurrence 
relation of hypergeometric function, 



d 
dz 



Fi (a; b;c;z) = ( — ) (a + 1; 6 + 1; c + 1; z) , 



where 



iFi {a,b;c;z) 



Tic) 



r{a + p)T{b + p) zP 

r(a)r(6)^ fx^T^) ^■ 



(53) 



(54) 



which is used to solve Eq. (22a) and present the corresponding lower component G„k(^) ^ 
follows 



1 



M + - a 

ro (M + - Cs) 



^riK 



d ^ _ ^ / X 
dr r 



-r/ro 



(K + H)ro 



+ 



' 1 - e-''/''o 

X 2F1 (-n, n + 2A„« + 27^^ + 1, 2A„« + 1; e"'^/'-") 

- (-n + 1, n + 2 (A„« + r?^ + 1) , 2 (A„« + 1) ; e"''/^") } . (55) 

Here, we remark that the hypergeometric series 2-F1 {—n, n + 2A„«; + 277^ + 1, 2A„k + 1; 
terminates for n = and thus converges for all values of real parameters t]^ and A„k- It is 
worthy to note that in the limit of ro goes to infinity then the Dirac spinor components 
become as 

lim Fr^^{r) = lim Gn«(r) = 0, 

ro— >-oo ro— >-oo 

which show that these components become unbound. 

For Cs > M+EnK and < M or < Af +EnK and E^k > M, we note that parameters 
given in Eq. (33c) turn to be imaginary, i.e., A^^ < in the s-state {k — —1) case. As a 
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result, the condition of existing bound states are Xnn > and 77^ > 0, that is to say, in the 
case of Cs > M + Enn and Enn < M, bound-states do not exist for some quantum number k 
such as the s-state {k = —1). Of course, if these conditions are satisfied for existing bound- 
states, the energy equation and wave functions are the same as these given in Eqs. (45), 
(52) and (55). 

In the spin-symmetric case, the numerical bound state energy spectrum is obtained from 
Eq. (40) with the choice of the following values for the parameters: M — 10 fnT^, tq — 10 
/m, Vq — 10 fmr^ and Cg — 10.1 fm~^. Thus the approximated values for the energy states 
using the usual approximation scheme (31) (i.e., with dg — 0) are presented in Table 2 for 
states n = 0,1,2 and 3 and spin-orbit quantum number k = ±1,±2,±3 and ±4 values. 
Apparantly, the spin unaligned states (k > 0), the orbital states (/ = |k, + 1/2| — 1/2, 
j = \k\ - 1/2): (lpi/2,0ci3/2), (2pi/2, lc?3/2, 0/5/2), (3pi/2, 2ci3/2, 1/5/2) and (34/2,2/5/2) are 
found to be degenerate in the presence of tensor (if 7^ 0) and in the limit of zero tensor 
couplings {H — 0). In addition, the energies of the npi/2, ndz/2 and n/5/2 states for 7^ 
are higher than the corresponding values for H — 0. Overmore, we reproduced the energy 
spectrum of Table 1 using the new approximation scheme (31) proposed in Ref. [44] through 
taking the dimcnsionlcss constant do ~ 1/12 [35,59-62]. The energy levels presented in Table 
2 are slightly different under our new approximation scheme. Further, the degeneracies in 
the above states have no Inoger exist. This is true for both Coulomb tensor {H ^ 0) and in 
the limit of zero tensor [H — 0) interactions. 

B. Pseudospin symmetric bound state solutions 

The exact pseudospin symmetry occurs in the Dirac equation when S{r) ~ —V{r\ In 
this part we will take the Hulthen potential for A(r), (i.e., = 0, or E(r) = Cyg — 

constant) , 



where Aq is a constant and the tensor potential as in Eq. (26). Thus, Eq. (24) can be 
reduced into the Schrodinger-type: 



A(r) 



(56) 




(e''Ao - 1) r\ 




(57) 
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where 

7 = Ajk (Afe - 1) , Ak^K + H, (58a) 

^^rl{En.-M-Cps)Vo, (58b) 

(K« + M) (K« - M - Cp,) < 7^0, En.<M + Cps, (58c) 

where k = —I and /t = / + 1 for /t < and k > 0, respectively. In the pseudospin symmetry, 
the eigenstates with with j — I ± ^ are degenerate for I ^ 0. The energy eigenvalues E^^ 
depend only on n and ie., = -E'(?T', '«('« — 1))- We now follow the previous procedures 
to obtain a differential equation for Gn^J^z\ 



{ {\-z) d 1 .2 1 1 ^ / ^ n 

\-d^-^W^)'dz-^^^^-^'''' -p,. + Po]jG„,(.) = 0, 



(59) 



where 

P2 = /5 + 7cio - ^n«, (60a) 

Pi=^+27do-7-2^nK, (60b) 

Po = l^Q - Snk = AL> (60c) 

and must be a positive real parameter for real solution. Thus, the energy spectrum is 
then given by 

. m^^±^. (61) 

2[n-\- K-\- H) 

The last formula gives an equation for the energy spectrum. When we insert Xnn and (5 into 
the above equation, then we find the following transcendental energy equation 

doiK + H){K + H - 1) 



(E^^ + M) {E^,-M-Cp,) 



'O 



rp (K« - M - Cp,) Fq _ (n + /€ + //) 
(n + K + i7) To 



2 



n = 0, 1,2, • • • ,n„,a. = roVH (£^n« -M-C,)-k-E, (62) 

where rimax being the largest integer which is less than ro-\/Vo (£^„« — M — Cs) — k — H, where 
> M -\- Cs- Obviously, the above equation can be reduced to Klein-Gordon solution 
when Cps — H — do — and after inserting q — 1, Vq — Sq and 5 — S± — I + 1, —I in Eqs. 
(8), (9) and (25) of Ref. [65]. Furthermore, the wave functions can be written as 

GnA^) = ^^^(l - z)^'^z^-- (-n, n + 2A„„ + 2Afe + 1, 2A„« + 1; -z) , 
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r(n + 2A„« + l) 

where 

, f 1+ + 1, K > 0, 
[ + K < 0. 

To avoid repetition, the solution of Eq. (24) can be found easily by applying appropriate 
paremetric transformations. A first inspection for the relationship between the present set of 
parameters {Snk, ^k) and the previous set {Snk, Vk) tells us that the energy solution, for 
pseudospin symmetry, can be obtained directly from those of the previous energy solution, 
in spin symmetric case, by applying the following appropriate parameter map [42,43]: 

FnA^) ^ Gn,,(r), ^ Cs ^ -Cps, V(r) ^ -V(r) (or Vq ^ -Vq), (65) 

or simply we apply the following transformations: 

FnAr) ^ Gn,n{r), M -M, Cs Cps, ^ A« (or k+H+1 k+H, k+H k+H-1), 

(66) 

on Eqs. (40) and (50) to obtain Eqs. (62) and (63) for energy spectrum formula and wave 
functions, respectively. 

IV. CONCLUDING REMARKS 

In this work, we have obtained analytically the approximate energy equation and the 
corresponding wavefunctions of the Dirac equation for the Hulthen potential coupled with 
a Coulombic-like tensor under the conditions of the spin and pseudospin symmetry using 
the parametric generalization of the NU method. The introduced Coulombic-hke tensor 
interaction generates additional centrifugal-like term L(L-|-l)r~^. Therefore, in order to solve 
the resulting Schrodinger-like equation analytically, a new approximation scheme was used 
to approximate the new spin-coupling centrifugal term 7r^^ arising from the Coulomb-like 
tensor which yields a closed form solution of the problem under consideration. The resulting 
solutions of the wavefunctions are written in terms of the orthogonal Jacobi polynomials or 
hypergeometric functions. Obviously, for exact spin when S{r) — V{r) {i.e., Cs — 0), 
the relativistic solution can be easily reduced to it's non-relativistic limit by the choice 
of appropriate mapping transformations. Also, in the limit of zero tensor couplings, the 

17 



present results reduce to the well-known solutions of the Dirac equation for the usual Hulthen 
potential with arbitrary spin-orbit coupling quantum number k. 
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Appendix A: Parametric Generalization of the NU Method 



We complement the theoretical parameterized formulation of the NU method for any 
arbitrary exponential potential by giving the essential polynomials, energy equation and 
wavefunctions together with their relevant constants as follows [42,43,63]. 

(i) The key NU polynomials: 

7r{z) = C4 + C5Z - [(v^ + C3v^) z - Vcsl , (Al) 

k=-{cr + 2C3C8) - 2v/^. (A2) 

r{z) = 1 - (c2 - 2C5) z-2 [(^ + C3^/c^) z-^/c^, (A3) 

r'{z) = -2c3 - 2 + C3v/^) < 0, (A4) 

(ii) The general energy equation: 

(C2 - C3) n + csn^ - (2n + 1) C5 + (2n + 1) + c^^) + C7 + 2c3Cs + 2^/^ = 0. (A5) 

(iii) The general wavefunctions: 

=z^i°(l-C3^r", (A6) 

cj){z) = z'^^'il - cszY'^ C12 > 0, ci3 > 0, (A7) 

yn{z) = P^<=iO'<=")(l - 2c3^), Mclo > -1, Men > -1, (A8) 

u{z) = ^fnZ^'^{l - C3zyPjf'°'''''\l - 2csz), (A9) 

where Pn'^\x), Mz/ > — 1, K/i > — 1 and a; G [—1, 1] are the Jacobi polynomials and J\fn is 
a normalization constants. Also, the above wave functions can be expressed in terms of the 
hypergeometric function as 

Unniz) = NnnZ"^^ (1 - Czzf^^ ^Fi {-U, 1 + Cw + Cn + Ti; Cio + 1; C3Z) , (AlO) 

where C12 > 0, C13 > and z e [0, 1/C3] . 

(iv) The relevant constants: 

1 1 

C4 = 2 (1 ~ ci) , C5 = - (C2 - 2C3) , cq = cI + Bi, 

C7 = 2C4C5 - B2, C8 = C4 + B3, Cg = C3 (C7 + CsCs) + Cq, 

2 

Cio = ci + 2c4 + 2^0^ - 1 > -1, Cll = 1 - ci - 2c4 H ^/c^ > -1, 

C3 

C12 = C4 + > 0, ci3 = -C4 H (v^ - C5) > 0. (All) 
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Appendix B: Normalization of the radial wavefunctions 



The normalization constant Muk can be easily calculated in closed form. To do this, we 
start by using the normalization condition [F„^(r)]^dr = 1, and under the coordinate 
change x — 1 — 2e~^/''°, the normalization constant Muk in (50) is given by 



n!r(2A„, + l) f fWl-xV^--' fl + x^^'"" 



T{n + 2A„« + 1) 



^p^2A„.,2^.)(^)j2^^^ (Bl) 



where A„«; > and rj^ > 1- The calculation of this integral can be done by writting one of 
the 2riK factors {1 + x) /2 in the form 



1 + x 



1 - 



1 — X 



2 2 ' 

and by making use of the following integral (see formula (7.391.5) in [67]): 

2^+T(n + iy+ l)r(n + + 1) 



(B2) 



J' {i-xy-^{i + xf [pi'''''\x)]^ dx ^ 



which is valid for Ru > and M/x > — 1. This leads to 



n\ur{n + u + + 



(B3) 



1 



r(2A„, + 1) 



2 An«r(n + 2A„« + l)r(n + 2A„« + +2ri^ + 1) 



To 



n\r{n + 2r]^ + 1) 



1/2 



(B4) 



20 



[1] J. N. Ginocchio, Relativistic harmonic oscillator with spin symmetry, Phys. Rev. C 69 (2004) 
034318-25. 

[2] J. N. Ginocchio, Pseudospin as a relativistic symmetry, Phys. Rev. Lett. 78 (3) (1997) 436-439. 
[3] J.N. Ginocchio, Relativistic symmetries in nuclei and hadrons, Phys. Rep. 414 (4-5) (2005) 
165-261. 

[4] P.R. Page, T. Goldman, J.N. Ginocchio, Relativistic symmetry suppresses quark spin-orbit 

splitting, Phys. Rev. Lett. 86 (2001) 204-207. 
[5] A. Arima, M. Harvey, K. Shimizu, Pseudo LS coupling and pseudo SU3 coupling schemes, 

Phys. Lett. B 30 (1969) 517-522. 
[6] K.T. Hccht, A. Adler, Generalized seniority for favored J 7^ pair in mixed configurations, 

Nucl. Phys. A 137 (1969) 129-143. 
[7] J.N. Ginocchio, D.G. Madland, Pseudospin symmetry and relativistic single-nucleon wave 

functions, Phys. Rev. C 57 (1998) 1167-1173. 
[8] A. Bohr, L Hamarnoto, B.R. Motelson, Pseudospin in Rotating Nuclear Potentials, Phys. Scr. 

26 (1982) 267-272. 

[9] J. Dudek, W. Nazarewicz, Z. Szymanski, G.A. Leander, Abundance and systematics of nuclear 
superdeformed states; relation to the pseudospin and pseudo-SU(3) symmetries, Phys. Rev. 
Lett. 59 (1987) 1405-1408. 

[10] D. Troltenier, C. Bahri, J. P. Draaycr, Generalized pseudo-SU(3) model and pairing, Nucl. 
Phys. A 586 (1995) 53-72. 

[11] J. Meng, K. Sugawara-Tanabe, S. Yamaji, P. Ring, A. Arima, Pseudospin symmetry in rela- 
tivistic mean field theory, Phys. Rev. C 58 (2) (1998) R628-631. 

[12] A.D. Alhaidari, H. Bahlouli, A. Al-Hasan, Dirac and Klein-Gordon equations with equal scalar 
and vector potentials, Phys. Lett. A 349 (2006) 87-97. 

[13] R. Lisboa, M. Malheiro, A.S. De Castro, P. Alberto, M. Fiolhais, Pseudospin symmetry and 
the relativistic harmonic oscillator, Phys. Rev. C 69 (2004) 024319-15. 

[14] H. Akgay, Dirac equation with scalar and vector quadratic potentials and Coulomb-like tensor 
potential, Phys. Lett. A 373 (2009) 616-620. 

21 



[15] H. Akgay, C. Tezcan, Exact solutions of the Dirac equation with Harmonic oscillator potential 
including a Coulomb-like tensor potential, Int. J. Mod. Phys. C 20 (6) (2009) 931-940. 

[16] M. Moshinsky, A. Szczepaniak, The Dirac oscillator, J. Phys. A: Math. Gen. 22 (1989) 
L817-L819. 

[17] V.I. Kukulin, G. Loyola, M. Moshinsky, A Dirac equation with an oscillator potential and 

spin-orbit coupling, Phys. Lett. A 158 (1991) 19-22. 
[18] G. Mao, Effect of tensor couplings in a relativistic Hartree approach for finite nuclei, Phys. 

Rev. C 67 (2003) 044318-12. 
[19] P. Alberto, R. Lisboa, M. Malheiro, A.S. de Castro, Tensor coupling and pseudospin symmetry 

in nuclei, Phys. Rev. C 71 (2005) 034313-7. 
[20] R.F. Furnstahl, J. J. Rusnak, B.D. Serot, The nuclear spin-orbit force in chiral effective field 

theories, Nucl. Phys. A 632 (1998) 607-623. 
[21] M.H. Moshinsky, Y. Smirnov, The Harmonic Oscillator in Modern Physics, Hardwood Aca- 
demic Publisher, Amsterdam, (1996) 289-404. 
[22] M.H. Pacheco, R.R. Landim, C.A.S. Almedia, One-dimensional Dirac oscillator in a thermal 

bath Phys. Lett. A 311(2003) 93-96. 
[23] J.Y. Gou, Z.Q. Sheng, Solution of the Dirac equation for the Woods-Saxon potential with 

spin and pseudospin symmetry, Phys. Lett. A 338 (2005) 90-96. 
[24] A. Arda, R. Sever, Approximate solution of the effective mass Klein-Gordon equation for the 

Hulthen potential with any angular momentum. Int. J. Theor. Phys. 48 (2009) 945-951. 
[25] Y. Xu, S. He, C.S. Jia, Approximate analytical solutions of the Dirac equation with the 

Poschl-Teller potential, J. Phys. A: Math. Theor. 41 (2008) 255302-255309. 
[26] Y. Zhang, Approximate analytical solutions of the Klein-Gordon equation with scalar and 

vector Eckart potentials, Phys. Scr. 78 (2008) 015006 (4pp). 
[27] W.C. Qiang, R.S. Zhou, Y. Gao, Application of the exact quantization rule to the relativistic 

solution of the rotational Morse potential with pseudospin symmetry, J. Phys. A: Math. Theor. 

40 (2007) 1677-1685. 

[28] A. Soylu, O. Bayrak, I. Boztosun, An approximate solution of Dirac-Hulthcn problem with 
pesudospiun and spin symmetry for any k state, J. Math. Phys. 48 (2007) 082302-9. 

[29] A. Soylu, O. Bayrak, I. Boztosun, k, state solutions of the Dirac equation for the Eckart 
potential with pseudospin and spin symmetry, J. Phys. A: Math. Theor. 41 (2008) 065308 

22 



(8pp). 

[30] C.S. Jia, P. Gao, X.L. Peng, Exact solution of the Dirac-Eckart problem with spin and pseu- 

dospin symmetry, J. Phys. A: Math. Gen. 39 (2006) 7737-7744. 
[31] L.H. Zhang, X.P. Li, C.S. Jia, Analytical approximation to the solution of the Dirac equation 

with the Eckart potential including the spin-orbit coupling term, Phys. Lett. A 372 (2008) 

2201-2207. 

[32] C.S. Jia, J.Y. Liu, L. He, L. Sun, Pseudospin symmetry in the relativistic empirical potential 

as a diatomic molecular model, Phys. Scr. 75 (2007) 388-393. 
[33] C.S. Jia, P. Gao, Y.F. Diao, L.Z. Yi, X.J. Xic, Solutions of Dirac equations with the Poschl- 

Teller potential, Eur. Phys. J. A 34 (2007) 41-48. 
[34] H. Motavali, Bound state solutions of the Dirac equation for the Scarf-type potential using 

Nikiforov-Uvarov method, Mod. Phys. Lett. A 24 (15) (2009) 1227-1236. 
[35] C.-S. Jia, T. Chen, L.-G. Cui, Approximate analytical solutions of the Dirac equation with 

the quantized Poschl- Teller potential including the pseudo-centrifugal term, Phys. Lett. A 373 

(18-19) (2009) 1621-1626. 
[36] X.-C. Zhang, Q.-W. Liu, C.-S. Jia, L.-Z. Wang, Bound states of the Dirac equation with vector 

and scalar Scarf-type potentials, Phys. Lett. A 340 (2005) 59-69. 
[37] X. Zou, L.-Z. Yi, C.-S. Jia, Bound states of the Dirac equation with vector and scalar Eckart 

potentials, Phys. Lett. A 346 (2005) 54-64. 
[38] G.F. Wei, S.H. Dong, Approximately analytical solutions of the Manning-Rosen potential with 

the spin-orbit coupling term and spin symmetry, Phys. Lett. A 373 (2008) 49-53. 
[39] K.-E. Thylwe, Amplitude-phase methods for analyzing the radial Dirac equation: calculation 

of scattering phase shifts, Phys. Scr. 77 (2008) 065005 (8pp). 
[40] A.D. Alhaidari, Solution of the Dirac equation by separation of variables in spherical co- 
ordinates for a large class of non-central electromagnetic potentials, Ann. Phys. 320 (2005) 

453-467. 

[41] C. Berkdemir, R. Sever, Pseudospin symmetry solution of the Dirac equation with an angle- 
dependent potential, J. Phys. A: Math. Theor. 41 (2008) 045302 (llpp). 

[42] S.M. Ikhdair, R. Sever, Approximate analytical solutions of the general Woods-saxon poten- 
tials including the spin-orbit coupling term and spin symmetry, to appear in Cent. Eur. J. 
Phys. (2010). 



[43] S.M. Ikhdair, Approximate solutions of the Dirac equation for the Rosen-Morse potential 

including the spin-orbit centrifugal term, to appear in J. Math. Phys. 51 (2) (2010). 
[44] S.M. Ikhdair, An improved approximation scheme for the centrifugal term and the Hulthen 

potential, Eur. Phys. J. A 39 (2009) 307-314. 
[45] S.M. Ikhdair, R. Sever, Any /-state improved quasi-exact analytical solutions of the spatially 

dependent mass Klein-Gordon equation for the scalar and vector Hulthen potentials, Phys. 

Scr. 79 (2009) 035002 (8 pages). 
[46] C.S. Lam, Y.P. Varshni, Energies of s eigenstates in a static screened Coulomb potential, 

Phys. Rev. A 4 (1971) 1875-1881. 
[47] B. Durand, L. Durand, Duality for heavy quark systems, Phys. Rev. D 23 (1981) 1092-1102. 
[48] R.L. Hall, Envelope representations for screened Coulomb potentials, Phys. Rev. A 32 (1985) 

14-18. 

[49] R. Dutt, K. Choudhury, Y.P. Varshni, An improved calculation for screened Coulomb poten- 
tials in Ray leigh-Schro dinger perturbation theory, J. Phys. A: Math. Gen. 18 (1985) 1379-1388. 

[50] J. Lindhard, P.G. Hansen, Atomic Effects in Low-Energy Beta Decay: The Case of Tritium, 
Phys. Rev. Lett. 57 (1986) 965-967. 

[51] I.S. Bitensky, V.K. Ferleger, LA. Wojciechowski, Distortion of H2 potentials by embedding 
into an electron gas at molecule scattering by a metal surface, Nucl. Instrum. Meth. B 125 
(1997) 201-206. 

[52] C.-S. Jia, J.Y. Wang, S. He, L.-T.Sun, Shape invariance and the supersymmetry WKB ap- 
proximation for a diatomic molecule, J. Phys. A: Math. Gen. 33 (2000) 6993-6998. 
[53] J. A. Olson, D.A. Micha, Transition operators for atomatom potentials: The Hilbert Schmidt 

expansion, J. Chem. Phys. 68 (1978) 4352-4356. 
[54] A.F. Nikiforov, V.B. Uvarov, Special Functions of Mathematical Physics, Birkhauser, Bassel, 
1988. 

[55] S.M. Ikhdair, R. Sever, Improved analytical approximation to arbitrary Z-state solutions of the 
Schrodinger equation for the hyperbolical potentials, Ann. Phys. (Berlin) 18 (10-11) (2009) 
747-758. 

[56] S.M. Ikhdair, Exact Klein-Gordon equation with spatially dependent masses for unequal 

scalar-vector Coulomb-like potentials, Eur. Phys. J. A 40 (2) (2009) 143-149. 
[57] J.D. Bjorken, S.D. Drell, Relativistic Quantum Mechanics, McGraw-Hill, NY, 1964. 

24 



[58] G.R. Satchler, Direct Nuclear Reactions, Oxford University Press, London, 1983. 

[59] Z.-Y. Chen, M. Li , C.-S. Jia, Approximate analytical solutions of the Schrodinger equation 

with the Manning-Rosen potential model. Mod. Phys. Lett. A 24 (23) (2009) 1863-1874. 
[60] Y.-F. Diao, L.-Z. Yi, T. Chen, C.-S. Jia, Arbitrary l-wavc bound state solutions of the 

Schrodinger equation with the Eckart potential. Mod. Phys. Lett. B 23 (2009) 2269-2279. 
[61] T. Chen, J.-Y. Liu, C.-S. Jia, Approximate analytical solutions of the Dirac-Manning-Rosen 

problem with the spin-symmetry and pseudo-spin symmetry, Phys. Scr. 79 (2009) 055002 

(7pp). 

[62] T. Chen, Y.-F. Diao, C.-S. Jia, Bound state solutions of the Klein-Gordon equation with the 

generalized Poschl-Teller potential, Phys. Scr. 79 (2009) 065014 (6pp). 
[63] S.M. Ikhdair, Rotation and vibration of diatomic molecule in the spatially-dependent mass 

Schrodinger equation with g-deformed Morse potential, Chem. Phys. 361 (1-2) (2009) 9-17. 
[64] S.M. Ikhdair, Bound states of the Klein-Gordon equation for vector and scalar general Hulthen- 

type potentials in L»-dimension, Int. J. Mod. Phys. C 20 (1) (2009) 25-45. 
[65] W.-C. Qiang, R.-S. Zhou, Y. Gao, Any Z-state solutions of the Klein-Gordon equation with 

the generalized Hulthen potential, Phys. Lett. A 371 (2007) 201-204. 
[66] F. Benamira, L. Guechi, A. Zouache, Comment on 'Exact solutions of the s-wave Schrdinger 

equation with Manning-Rosen potential', Phys. Scr. 80 (2009) 017001-3. 
[67] I.S. Gradshtein, I.M. Ryzhik, Tables of Integrals, Series and Products, Academic, New York, 

1969. 



25 



TABLE I: The NU parametric constants useful in calculating the energy eigenvalues and eigen- 
functions for the spin-symmetric Dirac equation. 



Constant Analytic value Constant Analytic value 



^1 


1 




Co 


1 


C3 


1 




C4 





C5 


-1/2 




C6 


92 + 1/4 


C7 


-qi 




C8 


\2 


eg 


{vk- 


1/2)' 


ClO 




Cll 




1 


Cl2 




Cl3 


Vk 
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TABLE II: The usual approximation to the relativistic energy spectrum (in fm^^) of the spin- 
symmetric Hulthen potential including the Coulomb coupling tensor of strength H = and H = 0.5 
for various values of n and k. 



n Z AC < En,^<o{H = or En,^^o{H = 0.5)« k>0 En,n>o{H = 0)» ^„,«>o(i? = 0.5)« 






-1 


OSl /9 






0. 


101446652 
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-1 


ISl /9 


0. 


,1057848200 


0. 


100388800 














2 


-1 


2Si /9 
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,1231107030 


0. 


136077055 
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-1 


351/2 
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,1518922615 


0. 


,170536938 
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-2 


0»Q/9 
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,101446652 
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Opi /9 


0. 


,123110703 


0. 


,1360770550 


1 1 


-2 
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,101446652 
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IPl /9 
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,1518922615 
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,1705369375 


2 1 


-2 


r Oj Z 
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,105784820 


0. 


100388800 


1 


ir V j A 


0, 


,191988313 


0. 


,2162202946 


3 1 


-2 
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0. 
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3jOl /9 
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,243203547 
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,2729055754 
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-3 
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2 


0/ z 


0. 


,1518922615 


0. 


,1705369375 


1 2 


-3 


14/2 


0. 


,105784820 


0. 


,101446652 


2 


14/2 


0. 


,191988313 


0. 


,2162202946 


2 2 


-3 


24/2 






0. 


,101446652 


2 


24/2 


0. 


,243203547 


0. 


,2729055754 


3 2 


-3 


34/2 


0. 


,105784820 


0. 


100388800 


2 


34/2 


0. 


,3052908168 


0. 


,3403207458 


3 


-4 


0/7/2 


0. 


,1518922615 


0. 


,136077055 


3 


0/5/2 


0. 


,191988313 


0. 


,2162202946 


1 3 


-4 


1/7/2 


0. 


,123110703 


0. 


,100388800 


3 


1/5/2 


0. 


,243203547 


0. 


,2729055754 


2 3 


-4 


2/7/2 


0. 


,105784820 


0. 


,101446652 


3 


2/5/2 


0. 


,3052908168 


0. 


,3403207458 


3 3 


-4 


3/7/2 






0. 


,101446652 


3 


3/5/2 


0. 


,3779539814 


0. 


,4181464025 



"We have taken rfo = during these calculations. 
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TABLE III: The new approximation to the relativistic energy spectrum (in fm^^) of the spin- 
symmetric Hulthen potential including the Coulomb coupling tensor of strength H = and H = 0.5 
for various values of n and k. 



n I 


K<0{l,j) 






> {1,3) En,n>o{H = or En,n>o{H = 0.5)« 





-1 


OSi /o 




0.1014318359 








1 


-1 


ISi /o 


0.105784820 


0.1004034771 








2 


-1 


2si /o 


0.123110703 


0.1360623873 








3 


-1 


3si /9 


0.1518922615 


0.1705222701 








1 


-2 


Od-j /o 
^x^o/z 


0.100057870 


0.1000005174 


1 


Opi/2 0.1232273855 


0.1362956215 


1 1 


-2 


-*-i^o/z 




0.1014893336 


1 


lpi/2 0.1014382946 


0.1707561995 


2 1 


-2 


2Dq /n 


0.100057870 


0.1130529719 


1 


2pi/2 0.1026480928 


0.2164398762 


3 1 


-2 




0.123227386 


0.1010361207 


1 


3pi/2 0.2433208010 


0.2731253525 


2 


-3 


Odr; /o 
0/Z 


0.123457217 


0.1132253472 


2 


04/2 0.1522417493 


0.17104639476 


1 2 


-3 




0.1061170006 


0.10163993177 


2 


W3/2 0.1923389117 


0.2167313327 


2 2 


-3 


24/2 




0.10163993177 


2 


24/2 0.2435547169 


0.2734174979 


3 2 


-3 


3(^5/2 


0.1061170006 


0.1132253472 


2 


34/2 0.3056423517 


0.3408332607 


3 


-4 


0/7/2 


0.1525865150 


0.13658291603 


3 


0/5/2 0.1926867981 


0.2171368804 


1 3 


-4 


1/7/2 


0.1000082780 


0.1135029417 


3 


1/5/2 0.2439041281 


0.2738248097 


2 3 


-4 


2/7/2 


0.0995403787 


0.09964377877 


3 


2/5/2 0.3059926539 


0.34124166206 


3 3 


-4 


3/7/2 




0.09964377877 


3 


3/5/2 0.3786566934 


0.4190685307 



"We have taken do = 0.0823058167837972 during these calculations. 
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